Phase randomness in a one-dimensional disordered absorbing 

medium 

V. Freilikher and M. PustilnikB 

The Jack and Pearl Resnick Institute of Advanced Technology, 
Department of Physics, Bar-Ilan University, Ramat-Gan 52900, Israel 

(February 1, 2008) 
Analytical study of the distribution of phase of the transmission coef- 
ficient through ID disordered absorbing system is presented. The phase is 
shown to obey approximately Gaussian distribution. An explicit expression 
for the variance is obtained, which shows that absorption suppresses the fluc- 
tuations of the phase. The applicability of the random phase approximation 
is discussed. 

PACS numbers: 42.25.Bs, 72.15.Rn, 05.40.+j, 78.20.Ci 



1 



Rigorous results for the reflection of waves from a disordered absorbing one-dimensional 
medium have been obtained a long time ago Nevertheless, despite its evident significance, 
the problem of wave propagation in absorbing disordered media has received less attention 
as compared to the studies of waves in passive (nonabsorptive) systems. In the pioneering 
work 0] the possibility of the localization of light has been predicted, which gave a new 
boost to the studies of classical waves in random media. It was believed during the years 
that the effect of absorption on the classical wave transport is similar to that of inelastic 
scattering processes on electrons in disordered solids. Indeed, it is well- known that when the 
incoherent effects are modelled by introducing an optical (imaginary) potential, the theory 
simplifies drastically 0. Only recently it was understood that the role of absorption differs 
fundamentally from the role of inelastic scattering for electrons |4j] (see also discussion in 
0, 0). The source of this difference is that absorption corresponds to the actual removal 
of particles, but does not affect the phase coherence. Conversely, the inelastic scattering 
causes the loss of the phase memory and suppresses the interference effects, but preserves 
the number of particles. Recently there have been a number of works dealing with the effect 
of absorption (and related problem of the coherent amplification) on the wave transport in 
random media ||- ]T2 . 



We study in this paper the probability distribution of the phase 9 of the transmission 
coefficient t through ID random absorbing media. It is well-established now that in a 
passive (nonabsorbing) ID medium the interference leads to the strong localization of all 
eigenstates in the presence of arbitrary weak disorder. If the length of the scattering region 
L is large compared to the characteristic length of elastic scattering /, the transmittance 
T = \t\ decreases exponentially with the length in average. On the other hand, T strongly 
fluctuates and obeys the log-normal distribution. If the disorder is weak, the phases of the 
reflection and transmission coefficients become random ||13|| , ||14|| . One may assume that 



the phase <p of the reflection coefficient is distributed uniformly in the interval (0, 2ir) and 
independently on the amplitude of the reflection coefficient and amplitude and phase of 
the transmission coefficient. This assumption will be referred below as the random phase 



approximation (RPA). RPA, however, does not implies that the phase 9 of the transmission 



coefficient is also uniformly distributed. Indeed, the recently obtained formula reads ||15|| : 

(t/t*) = (e 2ie ) = (l + L/l)e- 3L/l . (i) 



Note that calculating this average with the uniform distribution for 9 one would obtain zero. 

The arguments presented above suggest that RPA can be used for calculation of the 
distribution of the absolute phase 9 G (—00,00). When the absorption is zero, 9 relates to 
the density of states in the scatterer by p = n^ 1 [d(kL + 9)/dE], E = k 2 |n|. Therefore 



(k + 9/L) /tc is the total number of states per unit length, which is self-averaged in the 
limit L/l ^> 1 [[17]], and obeys almost Gaussian distribution with var(#/L) oc L" 1 . The 
natural question arises in the context of the discussion above is how absorption affects this 
distribution. This question is addressed below. 

A wave U(z) with the wavenumber k propagating in a one dimensional random medium 
satisfies the equation 

(J^ + k 2 e(zfjU(z) = 0, (2) 

where the dielectric permittivity e (z) is equal to 1 for z < and z > L and e(z) = 1 + 5e (z) 
for < z < L. We assume for simplicity that the imaginary part of Se(z), which is 
responsible for the absorption, is constant, while the real part e (z) is a random function of 
z, which may be approximated by Gaussian white noise: 

5t (z) =e{z)+ i/kla, (3) 

(£(*)> = 0, (s(z)s(z')) = ^5 (z-z'). (4) 

Here the length of elastic scattering / in a system without absorption and the length of 
absorption l a in a system without disorder have been introduced. The reflection (r) and 
transmission (t) coefficients are defined according to the wave behavior outside the disordered 
region: 




-ikz _|_ re %kz^ z > L 

z < 

To proceed, we introduce the joint distribution function 

V L (x, 6) = (5(x-x (L)) 5(9-9 (L))> , (5) 

where we have used the parametrization r = tanh (x/2) e l(p , t = \t\ e te , — oo < 9 < oo. To 
derive the Fokker-Planck equation for Vt (x, 9) it is convenient to use the so-called invariant 
imbedding equations for r and t (a detailed description of the similar calculation can be 
found for instance in M, Assuming that RPA is applicable for length L larger than 

the wavelength (kL ^> 1) and in the limit of weak scattering (kl ^> 1) and weak absorption 
(kl a ^> 1), we obtain the equation 

where 

Qp (x) = — In (sinhx) + /5coshx, (7) 



FW 4( 3 -^i?^) = i (2 + fls) - 



Here S = L/l and I is the scattering length for zero absorption. The parameter (3 = l/l a 
describes the relative strength of absorption. Alternatively, Eq.(O) can be derived by the 



method of Ref . |15| . Indeed, for j3 = Eq. (y) is equivalent to the equation for expectation 



values of Ref. [19|. It should be noticed that it is not necessary to assume the white noise 



distribution for e(z) to derive Eq.(^). It is enough to require the weakness of scattering 



kl > 1 |5|. 



After integration over 9 Eq. @ reduces to the well-known equation for the probability 
density distribution of x The solution of the later equation saturates at S ^> 4 to the 
limiting S'-independent distribution 0, || 

Poo (x) = /5e /3 - Q " (x) (8) 



with Q,p given by Eq. (|7|). The distribution (H) implies the following limiting value for the 
average reflectance 

f l-2/?m(l//3), /3<1 
i?oo= lim (i? 5 ) = l + 2/3e 2/3 Ei(-2/3) ~ <^ (9) 

( (2/T 1 , /? » 1 

X 

where Ei(x) = / dye y /y . 

—oo 

Here is a point to discuss the validity of RPA, which was essentiall for derivation of 
Eq.(||), and therefore underlies the result of Eq.fl^) that predicts a monotonic decrese of Rao 
with (3. Despite this result being known for more than twenty years, it was only recently 
understood ||, that this prediction becomes wrong for large (3. Moreover, it has been 
verified numerically that RPA fails at the same values of (3 ||. Let us estimate now the 
range of validity of RPA. In the limit (3 — l/l a ^> 1 the wave penetrates into disordered 
absorbing medium at a distance which is much smaller than elastic scattering length I. This 
means that disorder may lead only to a small correction to -Rqo, which can be taken into 
account by means of perturbation expansion in the small parameter 1/(3: 

r oc = r° + i^/r 1 + o (/r 2 ) . (io) 



The zeroth-order term here does not depend on (3 but on the parameter kl a only In the 
limit kl a <C 1 (strong absorption, which was the case in |J) this term is equal to 1. This is 
easy to understand since this limit corresponds to a large conductivity, where the skin-effect 
results in the almost perfect reflection. However, we are considering here the opposite limit 
of the weak absorption (kl a 1). In this limit, /3 > 1 corresponds to the absence of the 
elastic scattering rather than to the strong absorption, and R° is nothing but the reflectance 
of the ideal (ordered) lossy system. Calculating the reflectance Roo in the lowest order on 
(kla) -1 and we find that Eq. (TO) takes the form 

ik,~(4M a )- 2 +(2(3)- 1 [l + 0(l/kl)}. (11) 

Eq.flnp is exact for kl a ^> 1 and kl 3> 1. The first term in the r.h.s. in Eq.(|TT|) is the 
reflectance of the semi-infinite clean media with the wave number k\ = [k 2 (1 + i/kla)] 1 ^ 2 ~ 



k (1 + i/2kl a ) and is absent in Eq.(§). The second term represents the correction due to 
disorder and for (3 3> 1 coincides exactly with Eq.(^|). Since RPA was the only approximation 
made in derivation of the non-perturbative result Eq.(|9]), one immediately finds that RPA is 
applicable if, in addition to the weakness of both the scattering and absorption, the following 
inequality holds: (4kl a )~ 2 <C (2(3) _1 , or 

(3 < (3 ~ 8 (kl a f . (12) 

These simple arguments are consistent with the recent numerical simulations ||, which 
clearly show that RPA is valid only when (3 is smaller than some (3 (j3 <C (3^j, for which 
-Roo {(3) takes its minimum. Eq. ( |12"D provides therefore the value of (3 ~ (3$ in the leading 
order in the parameter 1/kl, characterizing the weakness of disorder. 

Now that we have established the limits of the validity of Eq.(|J) in Eq. (|T2|) , let us turn 
to its solution. It is natural to assume that correlations between the phase 9 and reflectance 
R are negligible, and distribution (El) factorizes (see e.g. p0|): 



V s (x,9)~Ps(x)P s (e), 5»1 (13) 



Integrating Eq. ([]) by making use of this conjecture, we get 



^ s(9) .i (2+(fls »| 



V s {ff) ^-(2+ (Rs)) -^Vs (ff) , (14) 



which immediately results in the Gaussian distribution with zero average (8) = and vari- 
ance 

var# = (1 + R OD /2)S. (15) 

One can see from Eqs.(p|),(|T5D that the absorption suppresses fluctuations of 6. Combining 
Eq.(|l5D with the recently obtained result var(lnT) = 2R 00 S |7j we obtain another useful 
formula 

var 9 = S + - var (In T) . (16) 
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Calculating (e 2%e j for (3 = with the Gaussian distribution specified by (15) one recovers 
Eq.(|l]) with an exponential accuracy. 

It is interesting to note that Eq.(fl5|) is rather general, and can be obtained without the 
conjecture ([HID. To do this we use the formula 



d 2 (e ipe 

/ p=0 

To calculate the characteristic function le tpe } we consider the Fourier-transform of 



J^M^te + ^H F K) Ps(i;p) ' (18) 

Integrated over x, P,s {x;p) gives us a characteristic function (e tpe ^. Following the procedure 
used in [0 we substitute 

P s (x;p) = e- n ? /2 ty(x;S) 

which yields an imaginary time Scrodinger equation for ty. Solution of this equation in 
spectral representation takes the form 



* (x; S 1 ) = £ e~ E " s ^ n (x) k (y) e n ^ 2 



where E n and ip n are eigenvalues and real normalized eigenfunctions of the Scrodinger op- 
erator 

( d 2 \ 

+ V p (x) ip n (x) = E n ip n (x) 



dx 



with potential 

1 p 2 (3 2 
V p (x) = — (l — sinh~ 2 H (3 — cosh -2 x/2^j — (3 cosh a; + — sinh 2 x (19) 



It can be easily seen [^] that the potential V p (x) has a minimum at x ~ ln(4//3), and 
the large 5 behavior of P$ (x; p) is dominated by the energy level localized in this well. 
For p = the energy of this level is exactly and the corresponding wave function is 
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ipQ = (Voo (x)) 1 ^ 2 = / 9 1 / 2 e /3 / 2 n p( x ) . For small p 2 (but for arbitrary (5) the perturbation 
theory yields: 



E (p) = {% P 2 F (x) ^) + o [p 2 ) = -p 2 (2 + R OQ )+o {p' 



Therefore for S ^> 4 



exp 



1 



-p 2 (2 + /U + o (; 



P 



Making use of Eq. (|17D , we end up with Eq. ([L5D . 

To conclude, we have studied analytically the influence of small absorption on the dis- 
tribution of the phases of the reflection and transmission coefficients through ID weakly 
disordered random media. It was shown that the phase if of the reflection coefficient can 
be treated as a uniformly distributed when the absorption is weak enough to satisfy the 
inequality Eq. fll~2]) . An explicit expression for the variance of the phase 9 of the transmission 
coefficient was obtained (Eq.flToD). The result shows that the presence of the absorption 
slows down the evolution of 9 towards a random limit. 

We thank I. Yurkevich for helpful discussions. 
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